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ABSTRACT 

The  abundance  anomalies  associated  with  lie-rich  flares,  are 
attributed  to  the  Intrinsically  anomalous  acceleration  behavior  of 
species  with  ^  <  1,  ,  in  the  presence  of  hydrogen  cyclotron  waves. 
IThile  the  overall  model  Is  along  the  lines  of  Fisk's  (1978)  proposal, 
the  difficulties  associated  with  triggering  of  nonhydrogenlc  cyclotron 
waves  are  avoided. 
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l.  INTRODUCTION 

One  of  the  most  intriguing  results  of  energetic  sol.tr  particle  flux 
measurements  at  ~  1  AM  Is  the  discovery  of  anomalous  enhancements  In  the 
abundance  of  some  ionic,  species  during  occasional  "unusual"  events, 
often  called  'lie-rich  flares  (Ramaty  et  al.  1980  and  references 
therein).  Three  different  types  of  enhancements  can  he  distinguished  in 
these  events: 

O 

t)  the  enhancement  of  the  isotope  'lie  hy  ns  much  as  3.3  orders 
of  magnitude,  at  energies  —  1-10  Me V/ nucleon. 

II)  enhancement  of  heavy  (A  >  '< )  tons  by  a  factor  of  5-10. 

Hi)  Deviations  of  other  charge  states  from  solar  abundances. 
These  are  often  correlated  during  ^lie-rich  flares.  The  most  promising 
theoretical  models  consider  a  two  stage  acceleration  process  (Ibragimov 
et  al .  1973,  Fisk  1973);  a  preacceleration,  or  preheating  stage  due  to 
wave-particle  Interactions,  followed  by  a  process  requiring  a  threshold 
injection  velocity  such  as  a  Fermi  type  acceleration  (Mbbius  et  al. 
19R2). 

Of  the  above  models  the  one  proposed  by  Fisk  (1973)  seems  to  be  the 
most  successful  in  explaining  not  only  enhancement  type  i)  but  types 

a 

it)  and  ill)  as  well.  It  is  based  on  selective  preheating  of  the  lie 
and  certain  heavy  Ions  such  as  Fe  by  resonant  interaction  with  ion 
cyclotron  waves.  Such  waves  can  be  driven  unstable  hy  electron  currents 
or  Ion  beams.  In  a  pure  hydrogen  plasma  the  ton  cyclotron  waves  have 

O 

frequency  above  the  proton  cyclotron  frequency  (l.e.  u>  »  (,*“  + 

2  1  1/2 

k  c“  )  =  1 . 2—  1 . A >*., .  c„  Is  the  sound  speed)  (Kadomtsev  1965), 

Is  11  n 

However,  In  the  presence  of  notlcable  amount  of  doubly  Ionized  and 

for  some  combination  of  drift  velocities  and  plasma  parameters,  waves 


f  I  , 

np.ir  the  cyclotron  frequency  can  be  excited.  Since 

a)  —  l .  2i2  =  12  ,  ^llo++  can  be  resonantly  accelerated.  In 

‘4lle  +  *  A 

addition,  ions  with  —  -  3.3  and  —  ~  4.5  can  be  accelerated  by  the 

harmonics  (i.e.  ui  =  ni-j“k.2v2)  although  less  effeciently. 

The  subject  of  ion  energization  by  ion  cyclotron  waves  in 

multispecie  plasmas  has  been  recently  examined  by  Papadopoulos  et  al . 

(1980),  Singh  et  al.  (1981,  1982),  and  Varvoglis  and  Papadopoulos 

(1982).  It  was  noted  that  for  wave  amplitudes  exceeding  a  threshold, 

* 

which  depends  on  — ,  ion  acceleration  does  not  require  resonance,  but 
can  be  very  strong  even  for  significant  frequency  mismatch.  The 
physical  process  responsible  for  this,  is  resonance  overlapping.  The 
general  theory  was  given  by  Chirikov  (1979),  while  its  application  to 
lower  hybrid  heating  of  fusion  plasmas  in  several  publications  (Fukuyama 
et  al.  1978,  Karney  and  Bers  1977,  Hsu  1982).  It  has  been  called  non¬ 
resonant  stochastic  acceleration  to  be  distinguished  from  the  one 
requiring  cyclotron  resonance  (i.e.  u>  =  ni2^  -  k.,vz))  and  as  first  shown 
by  Papadopoulos  et  al.  (1980)  exhibits  a  large  selectivity  in  the  ionic 
A/Q. 


In  this  paper  we  revise  the  model  advanced  by  Fisk  (1978)  by 
including  the  proper  non-1 inear  physics  of  particle  energization  by 
electrostatic  ion  cvcltron  waves  (EIC).  Our  model  retains  two  basic 
concepts  of  Fisk  (1978).  The  energization  by  EIC  waves  and  the  need  for 
a  second  stage  acceleration.  However  there  is  no  need  for  exciting 

A  j  t- 

He  cyclotron  waves,  since  the  dominant  process  is  non-resonant  and 

A 

can  be  accomplished  by  hydrogen  cyclotron  waves.  The  —  selectivity  in 
the  flux  available  for  energization  in  the  second  stage  process,  enters 

fit 

through  the  non-linear  saturation  level  ~L 

TH 


,  which  In  conventional 
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theories  ( Out')  and  Dupree  107'),  Paltnadesso  et  al,  \cn depends  on  the 
current  that  drives  the  instability. 


1 1  *  N0N-LINE\R  STOCHASTIC  ACCKI.FKATION 

The  non-linear  energization  of  minority  species,  taken  as  test 
particles,  by  an  F.IC  wave  with  typical  frequency  <o  near  and 

propagating  perpendicular  to  the  ambient  magnetic  field  can  be  described 
as  a  diffusion  in  velocity  space,  even  in  the  absence  of  resonance  (i.e. 
m-k^v  t  n  ’ji)  (Ahe  et  al.  (1981),  Hsu  (1982),  Varvoglis  and  Papadopoulos 
(1982)),  if 

’  R(v)  r  (1) 


ei-  s  1_ 

T  -  U  2 

K  <fc,y  »->„(«•) 


In  eqn.  (1)  >;>  is  the  wave  potential,  v  =  --  n  =  |  +  .‘5]  (the 

v,  'll  ’ 

integer  nearest  to  v) ,  r  =  (k^R^)  ^  ~,  Jn  is  the  ncl)  order  Bessel 


function  of  the  first  kind,  ind  R^,  T,j,  the  gvroradius,  temperature 
and  thermal  velocity  of  hydrogen,  which  is  assumed  to  bo  the  dominant 
species.  R( v)  is  a  function  whose  details  have  been  determined 
computationally  by  Hsu  (1982),  which  for  v  <  3  oscillates  around  R=1  and 
approaches  R=l  for  >  3.  (Note  that  vR(v)  =  P(  >)  is  plotted  in  fig.  5 
of  Hsu  (1982)1.  From  (1),  we  can  find  the  upper  and  lower  limit  of 

O  V 

stochast icitv  In  velocity  space  for  a  given  —  value  as  done  In 

Tn 

Papadopoulos  et  al.  (19801.  Here  we  take  a  different  approach  and  vise 

C  t1 

(l)  to  determine  the  minimum  value  of  (■=—)  required  for  the  lower 

thr 

velocity  stochastic ity  limit  to  be  equal  or  below  the  thermal  speed  Vj 

C  Z  t*  t> 

of  a  particular  specie  j.  This  impltes  that  for  =—  >  (,=— )  the  bull 

11  Si  thr 

of  the  species  will  be  preheated.  Notice  that  for  an  isothermal  plasm, 
(by  Isothermal  here  we  refer  only  to  ion  temperatures) 


•Maui 
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r 


.1/2  V, 

(kiV  -o-  V- 
) 


(2) 


T!ie  value  of 
equality  and 


(— — )  ran  thou  ho  computed  from  (1)  by  considering  the 

It  thr  V 

solving  for  ~  *  I,  which  corresponds  to 

j 


r 


r 

J 


<k,'V 


O 


(3) 


(A  and  O  correspond  to  specie  j).  To  he  concrete  we  examine  F1C  waves 


with  to  -  1.2  and 


(k  R^)‘  =  1.5  ,  wiildi  correspond  to  maximum  linear 


growth  (Kadomtsev  ,  1965).  it  is  easy  to  see  that  for 

>">  I  which  is  of  course 
unrealistic.  We  can  safely  conclude  that  for  A/Q  >  2.5  ,  only  tail  or 


•\ 

v  >  3.5,  which  corresponds  to  -jr  >  2.5,  (— -) 

)!  thr 


supratherraal  particles  can  be  accelerated  in  an  isothermal  plasma. 
However  in  the  region  of  small  —  the  situation  changes.  This  is  seen 


C*  I1 

from  fig.  I,  where  (— ) 

A  '*  tlir 

v  *  —  —  for  v  <  3.  Two  important  effects  are  shown  la  fig.  1.  First, 
**H  Q 

In  the  region  of  =  1.5-2,  which  Involves  ^!le++  and  ~tlle++, 

G  i* 

(~ — )  <  1,  so  that  bulk  preheating  of  these  species  is  possible. 

II  thr  s 

Second  the  threshold  is  a  non-mono tonic  function  of  — —  and  —  (e.g. 

"  H  • 

for  —  =  1.2,  (.p-)  =  .40  for  £=1.3  and  .45  for  £=  2.0). 

'll  Ml  thr  }  1 

These  results  were  confirmed  In  a  tost  particle  simulation  of  100 


Is  shown  as  a  function  of 


tons  with  an  initial  Maxwellian  distribution  with  T^  =  Tjj  interacting 
with  hydrogen  cyclotron  waves  at  w  *  1.2  W..  with  amplitude  —  =  0.45 


H 


propagating  perpendicularly  to  a  homogeneous  static  magnetic  field. 

This  amplitude  Is  consistent  with  auroral  observations,  laboratory 
experiments  ( Ktutner  et  al.  1973,  Ulihmer  1976)  as  well  as  theoretical 
estimates  (Dum  and  Dupree  1970,  Palmadesso  et  al .  1974).  The  result  for 


7 


^lle^  ions  is  shown  in  Fig.  2  for  t  =  0  and  t  ~  170  .i  '  .  Note  that  in 
tho  case  shown  in  Fig.  2  a  value  Ak  =  -—  lias  boon  used  as  the  wnve 
spread . 

From  Fig.  2  we  see  that 

i)  apparently  the  lover  velocty  threshold  rm^n  is  in  this 

case  nearly  zero,  since  the  whole  distribution  seems  to 
Interact  with  the  wave 

111  the  maximum  velocity  is  r  *  7  which  corresponds  to 

max 

-  7  v_  “4  v„,  where  vn  is  the  thermal  velocity  of 
3lle  ” 

3,do 

iil)  the  acceleration  time  scale  is  t  ■*  150  .!  '  ~  0.2  msec 

o  4*4-  „ 

iv)  a  considerable  number  of  JHe  ions  (~  95%)  have 
velocities  greather  than  v-, 

tie 

Note  tliat  the  energy  AK  absorbed  by  tlie  lie  ions  in  this  acceleration 

process  is  trivial;  e.g.  for  interaction  region  of  the  order  of  10  kn  , 
l7  13 

AE  ~  10  “  if  ergs. 

The  same  parameters  were  used  in  a  run  with  Sle^  ions.  Their 
distribution  function  was  unaffected  at  all  times,  e.onsistent  with  the 
fact  that  the  wave  amplitude  was  below  threshold  for  this  species. 

t  4*4- 

According  to  our  scenario,  almost  all  lie  ions  (but  only  the  Ions  in 
the  tall  of  the  ^llo+t  distribution)  will  be  injected  to  the  second  stage 
acceleration.  The  ^He+^/"*lle++  ratio  of  the  injected  Ions  will  depend  on 
the  exact  profile  of  tho  second  stage  threshold  function  (e.g.  its 
steepness),  but,  as  shown  by  Fisk.  (1978),  this  ratio  can  easily  attain 
values  of  order  unity.  if  tor  Instance  the  threshold  function  is  a 
steep  function  it  /  *  3v  then  ^  80%  of  the  He  ions  In  Fig.  2  wil7 

be  Injected,  In  contrast  with  only  0.0002%  of  the  He  ions.  Moreover, 
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If  the  space  averaged  wave  amplitude  fluctuates  with  time,  intermediate 
values  of  the  *lle^+/Mlle t-f  ratio  could  result,  from  l  (amplitude 
consistently  between  0.40  and  0.45)  to  the  normal  ~  10  4  (amplitude 
consistently  above  .45).  Note  that,  because  the  time  scale  for  this 
acceleration  mechanism  is  fast,  even  irregular  msec  fluctuations  of  the 
wave  amplitude  above  threshold  can  efficiently  contribute  to  particle 
acceleration,  while  grad  B  and  finite  kz  effects  are  not  expected  to 
affect  tlie  results. 


IV.  HEAVY  IONS 

The  previous  section  focused  on  the  He/  He  anomalies  during  He- 
rlch  events.  In  this  section  we  discuss  the  implications  of  the 
mechanism  to  the  acceleration  of  other  ionic  species;  specifically  the 


anomalies  on  various  charge  states  of  C,  0  and  Fe  ions,  associated  with 

^He-rich  events  (e.g.  Anglin  et  al . ,  1977;  Mason  et  al.,  1979;  Ma  Sung 
et  al.,  1981). 

a)  C  ions.  The  depletion  of  C  in  some  lie-rich  events  (Mason  et 
al.  1979)  can  be  accounted  by  the  same  process  that  creates  the  Z,He++ 


depletions  discussed  before,  since  for  the  typical  corona  temperatures 
of  ~  2  x  106  °K  almost  all  (>  9(H)  of  C  ions  are  in  the  +6  charge  state, 


having  the  same,  as  the  4!lo"f~1’  ions. 


A=  , 


For  the  rest  of  the  C  ions 


the  wave  amplitude  threshold  does  not  necessarily  increase,  so  that  C 


Lons  whith  Q  /=  6  can  be  non-resonnntly  accelerated.  As  a  result  while 
the  overall  C  flux  In  C-poor  flares  is  low  (~  of  the  normal,  If  we 

use  the  relative  abundance  of  C+<4  as  a  guideline)  almost  all  of  it 
+5 

should  be  In  the  C  state  since  this  is  the  next  more  abundant  after 
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b)  Fe  Lons:  Mci  Sting  et  .il  .  (1981)  and  Kleoker  et  si l  .  (  1982)  found 

that  In  at  least  one  Sle-rich  event  Ke  was  in  a  higher  than  usual  mean 

charge  state  (i.e.  Q  =  I'O,  with  Ions  having  sis  large  Q  as  13.  These 
\ 

ions  have  ■—  >  3  and  no  major  flux  anomalies  would  he  expected  if  their 
temperature  is  equal  to  Tl( .  If,  however,  Llieir  temperature  Tp  >  6TU , 


values  of  (— — )  .4-. 3,  will  result  in  accelerating  a  considerable 

T11 

fraction  (2,  10%)  of  the  available  states  o£  Fe  ions,  as  can  be  seen  from 

eqn.  (1).  If  we  assume  that  T  »  10**  UK,  large  enhancement  in  Fe , 

ions  would  require  T  N  b  x  10**  in  which  case  more  than  99%  of  Fe 

le 

will  be  in  charged  states  with  Q  >  IS.  In  our  model  we  associate  the 

enhancement  of  Fe  1  ,  with  the  presence  in  the  first  acceleration  region 

of  Fe  with  T  >  hi  ,  potentially  converted  from  hotter  regions. 

Fe  —  II 

c)  0  ions:  For  these  ions  also  >  3,  and  the  lower  velocity 
threshold  for  acceleration  will  he  the  same  for  0+1  and  0+f’ .  However 
the  maximum  acceleration  energy  (Papadopouios  et  al ,  1980)  scales  as 

Q"-^~,  so  that  larger  Q+:>  fluxes  can  enter  the  second  acceleration 

stage  than  for  the  same  initial  states.  This  could  cause  an 

+S  1 

occasional  0  enhancement  in  He  rich  flares  as  seen  by  ILa  Sung  et  al. 
(1981). 


V.  DISCUSSION  AND  CONCLUSIONS 

In  tills  paper  we  tiave  attempted  to  attribute  the  abundance 
a 

anomalies  observed  in  lie-rich  events,  to  the  intrinsically  anomalous 

acceleration  behavior  of  minority  species  with  ^  <  3,  in  the  presence 

of  hydrogen  cyclotron  waves.  The  theory  was  based  on  the  observations 

—  <  3 

Q 


that  for 
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(1)  Tin*  threshold  for  buik  non-resonant  acceleration  is  relatively 

low  (i.o.  — ' -  .A-.  5),  and  it  is  therefore  possible  to  have  Sie/JI  and 

TH 

a  , 

lie/H  enhancements. 

el 

(2)  Snail  variations  in  — —  or  the  wave  frequency  w,  can  produce 

rH 

major  variations  in  the  accelerated  fluxes  resulting  in  erratic  behavior 
in  the  abundance  ratios  of  ions  with  —  <  1. 

'I 

The  present  theory  improves  on  the  model  presented  by  Fisk  (197R) 
by  avoiding  the  difficulties  associated  with  the  special  conditions 
required  to  excite  ~‘lle++  cyclotron  waves  and  its  harmonics  at  large 
amplitudes,  and  treating  properly  the  non-linearity  associated  with 
large  amplitude  acceleration.  We  should  caution  the  reader  that  we  have 
examined  only  one  aspect  in  the  chain  of  the  acceleration  events.  A 
complete  picture  requires  input  on  the  plasma  composition  and 
temperatures  at  the  preacceleration  site  and  a  description  of  the  second 
stage  acceleration  processes.  In  this  respect  any  observational  support 
of  our  theory  might  be  clrcumstanctal .  We  feel,  however,  that  the 
coincidence  of  the  observed  ^  region  of  anomalies,  with  the  region  of 
anomalies  expected  on  the  basis  on  intrinsic  stochasticlty  behavior  for 
hydrogen  cyclotron  waves  cannot  bo  coincidental. 
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